The non-integrability of the Zipoy-Voorhees metric 
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The low frequency gravitational wave detectors like eLISA/NGO will give us the opportunity 
to test whether the supermassive compact objects lying at the centers of galaxies are indeed Kerr 
black holes. A way to do such a test is to compare the gravitational wave signals with templates of 
perturbed black hole spacetimes, the so-called bumpy black hole spacetimes. The Zipoy-Voorhees 
(ZV) spacetime (known also as the 7 spacetime) can be included in the bumpy black hole family, 
because it can be considered as a perturbation of the Schwarzschild spacetime background. Several 
authors have suggested that the ZV metric corresponds to an integrable system. Contrary to this 
integrability conjecture, in the present article it is shown by numerical examples that in general ZV 
belongs to the family of non-integrable systems. 
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1. INTRODUCTION 

It is possible that in the next decade the European 
New Gravitational Wave Observatory (NGO) P, [3| will 
be launched. The NGO will allow the tracing of the 
spacetime around supermassive compact objects (10^ — 
10^ Mq ), lying at the centers of galaxies, by detecting the 
gravitational waves emitted by an inspiraling much less 
massive compact object, e.g. a stellar mass black hole 
or a neutron star. The motion in such a binary system 
is known as Extreme Mass Ratio Inspiral (EMRI). We 
anticipate that these supermassive compact objects are 
black holes and the spacetime around them is described 
by the Kerr metric. 

However, we must check this hypothesis. An experi- 
mental test is to extract information from gravitational 
waves emitted by an EMRI system. Ryan showed [3] 
that we can extract the multipole moments of the cen- 
tral body from the gravitational wave signal (see also 
0, IE])- Thus, any non-Kerr multipole moments should 
be encoded in the waves. This was basically the idea 
Collins and Hughes originally conceived and stated in 
0. They, therefore, constructed a perturbed black hole 
which remained stationary and axisymmetric, and named 
it a "bumpy" black hole. According to their original idea, 
we should be able to measure deviations of the multipole 
moments of this perturbed black hole from those of a Kerr 
black hole in cases of EMRIs. But instead of following 
this approach, they preferred to explore the bumpiness 
of the perturbed spacetimes by measuring the periapsis 
precession Q . Their original idea was implemented later 
by Vigeland Q- 

Since Collins and Hughes' work, various approaches 
[isj , have been applied in order to identify the imprints of 
a perturbed spacetime background on the gravitational 
wave signal in the case of an EMRI system. On the 
other hand, the gravitational wave spectra are not the 



'Electronic address: ' gglukes®gmail.com| 



only spectra which can probe the bumpiness of a black 
hole. Electromagnetic spectra have been suggested in 
several works p6 19] as a plausible alternative. For more 
information about bumpy black hole detection methods 
see [13, mi- 

In previous works, we approached the bumpy black 
hole topic by studying the non-linear dynamics of the 
geodesic motion on the corresponding perturbed back- 
ground m, [12, H^l . Phenomena related to non-linear dy- 
namics appear in the bumpy black hole spacetimes due to 
the fact that the corresponding (stationary and axisym- 
metric) systems are missing a fourth integral of motion 
analogous to that of the Kerr spacetime, i.e. the Carter 
constant [13]. In fact, the existence of this extra sym- 
metry, associated with the Carter constant, is the reason 
that the Kerr system is integrable. 

The Carter constant seems to hold only for certain 
types of systems (see [2^, [2^ and references therein) 
which include the Kerr system. However, it is not typical 
to find spacetimes which are integrable. 

The Zipoy-Voorhees (ZV) metric [13,[2l|, known also 
as the 7 metric |20|, has been conjectured to correspond 
to an integrable system. In particular, in [s^l the authors 
examined the ZV spacetime numerically and did not find 
signs of chaos. A new numerical investigation by Brink in 
[31| gave the same result, which led Brink to conjecture 
that the ZV spacetime is integrable. Hence, she tried to 
find the missing integral of motion |3l| . The question 
of the "missing" integral was recently addressed also in 
(3^ . In this article the authors initially conjectured that 
the ZV is integrable in order to search for the missing 
integral, however their investigation led them to the op- 
posite direction. Namely, they proved the nonexistence 
of certain types of integrals of motion in a specific case 
of the ZV metric. 

The integrability conjecture is certainly true in two 
special cases, i.e. when a "free" parameter of the metric 
gives the Minkowski metric or the Schwarzschild metric. 
But, in the present article it is shown that this does not 
hold in general. In general, the phase space of the ZV 
system has all the features of a perturbed system, like 
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chaotic layers, Birkhoff chains etc. These features are 
found by using Poincare sections and by employing the 
so-called rotation number indicator |33l - l35l |. Thus, by 
numerical examples it will be shown that in general the 
ZV metric corresponds to a non-integrable system. 

The paper is organized as follows. Section [2] summa- 
rizes some basic elements regarding the ZV metric, the 
geodesic motion in the ZV spacetime background, as well 
as some essential theoretical elements regarding Hamil- 
tonian nonlinear dynamics. Numerical examples of the 
ZV non-integrability are presented in section[21 Section^ 
summarizes and discusses the main results of the present 
work. The accuracy of the integration method used to 
calculate the geodesic orbits is discussed in appendix El 



2. THEORETICAL ELEMENTS 

2.1. The Zipoy-Voorhees spacetime 

The Zipoy-Voorhees metric [23, describes a two pa- 
rameter family of Weyl class spacetimes, which are static, 
axisymmetric and asymptotically flat vacuum solutions 
of the Einstein equations. The ZV line element in pro- 
late spheroidal coordinates is given by 

ds^ = gtt dt^ + dx^ + gyy dy^ + g^^ dcj)^ , (1) 

where the metric elements are 

gtt 



where 
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and k — M/5 is the ratio of the source mass M and 
of the "oblateness" parameter 5. The ZV quadrapole 

moment is Q = Af^^^^— ^ [s^. If 5 = 1 the ZV metric 
3 

([2]) describes the standard Schwarzschild spacetime and 
the corresponding central object is spherically symmetric 
(2 = 0). If (5 > 1 the ZV metric describes a spacetime 
around a central object more oblate than a Schwarzschild 
black hole; if < 5 < 1 the central object is more prolate. 
Finally, if (5 = we get the Minkowski flat spacetime. 

The ZV metric takes a more familiar form [2^[4^ by 
using the transformation 



y = cos 9. (3) 

Then, the metric element ^ is 

ds^ = -Fdt^ + F-^[G dr^ + H d6^ 

+ (r^ - 2k r)siT?9d(t?], (4) 



F = 



G 



H = 



2k r 



-2k r + k^ sin^ 9 
(r2 - 2k rf 



(5) 



From this formulation, it is easy to check that for S = 
1 we get the Schwarzschild spacetime around a central 
object of mass k This means that the parameter S 
also "measures" how much more (or less) mass M = S k 
the ZV central object has compared with a Schwarzschild 
black hole of mass M — k. 

For < (5 ^ 1 the event horizon is broken and 
a curvature singularity appears along the line segment 
p ~ 0, |z| < fc [37]. Thus, the ZV spacetimes in general 
describe naked singularities (for further information see 
[3§-38]). On the other hand, there are no signs of closed 
timelike curves {g^^ > 0, V x, y). 

The parameter S has one more property of present in- 
terest: from a dynamical point of view S can be seen 
as a perturbation parameter of the Schwarzschild space- 
time. When the value of 6 departs from 1, the spherical 
symmetry is broken and the new axisymmetric spacetime 
needs not possess a Carter-like constant. In fact, some 
indication of this absence was found in 32], where for 
the 5 — 2 case the authors proved that a broad category 
of integrals of motion has to be excluded from being the 
integral that would make ZV integrable. In section [3] it is 
shown that indeed such an integral cannot exist, because 
ioi S — 2 there are geodesic orbits which are chaotic. 
Similar results are found for other representative values 
of 5 as well. 



2.2. Geodesic motion in the Zipoy Voorhees 
spacetime background 

The equations of geodesic motion of a "test" particle 
of rest mass /i in a spacetime given by the metric g^i, are 
produced by the Lagrangian function 



L 
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(6) 



where the dot denotes derivation with respect to proper 
time and the Greek indexes stand for spacetime coordi- 
nates. The Lagrangian ^ has a constant value L — 
along a geodesic orbit, due to the g^i, i^i" = — 1 con- 
straint. 

Since the ZV spacetime is axisymmetric and stationary 
the corresponding momenta 



= dL/dx" 

are conserved. These are the specific energy [48 
dL 



E 



di 



/m = -gttt, 
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FIG. 1; (a) The black curve shows the effective potential VcB for 5 = 2, k — 1, E = 0.97, Lz = 7.6, z — along the p axis. 
The letters A, C, D indicate the positions of the Vca roots for z = 0, while the black horizontal line shows the zero value of 
the VcS axis. The letter B indicates the position of a local maximum of the Ves, which separates the plunging orbits from the 
bounded non-plunging, (b) The effective potential cut along the VcB — plane. The A, C, D show the same as in (a) and the 
CZV corresponds to the curves of zero velocity. 



and the specific azimuthal component of the angular mo- 
mentum 

dL 

-^-z = -^/a* = 500'/' ■ (9) 

For brevity, these two integrals are hereafter simply re- 
ferred to as the energy E and the angular momentum 

The z index refers to the cyhndrical (p, z) coordinate 
system used in all figures. The transformation equations 
from spheroidal prolate to cylindrical coordinates is 

p = fcv/0^^~^I)O~y^, z^kxy . (10) 

The set of coordinates (p, z) describes the meridian plane. 
Due to the two integrals of motion (|5]), © we can restrict 
our study to this plane. We just have to reexpress ([8]), 
(O to get i and 4> as functions of E and L^, i.e. 

i = -E/gtt, = L^/g^^ , (11) 

and then substitute t, (j) into the two remaining equations 
of motion. Then, from the original set of 4 coupled second 
order ordinary differential equations (ODEs), we arrive 
at a set of 2 coupled ODEs. 

The motion on the meridian plane satisfies yet another 
constraint. If we use the transformation ([TU|) on the met- 
ric (HI) , substitute ([TT]) and replace the Lagrangian func- 
tion with its constant value, we get 

l{p^ + z^) + V,s{p,z)^0, (12) 



where 

1 / E'^ \ 

V,s - -Z 1 + — -t- ^ , (13) 

2 V 9tt 94,4, J 

r± = ^p^ + {z±ky. 

The VcS corresponds to a Newtonian-like two- 
dimensional effective potential (Fig. [1]). The roots 
{VcS — 0) of this effective potential (Fig. [1^) produce 
a curve on the meridian plane (p, z) called the curve 
of zero velocity (CZV) (Fig. [TJ)). This nomenclature 
follows from noting that whenever Vcs — the constraint 
()12p gives p = i = 0, which means that whenever a 
geodesic orbit reaches the CZV the velocity component 
in the p, z plane becomes equal to zero. 

On the meridian plane the CZV determines the set 
of initial conditions of geodesic orbits not escaping to 
infinity. These bounded orbits can either plunge to the 
central ZV compact object or revolve around it. In Fig. 
[1^ the plunging orbits lie before point A, while the non- 
plunging orbits lie between the points C and D. In Fig. [T]d 
the point A lies on the CZV containing all the plunging 
orbits, while C and D lie on the CZV containing the non- 
plunging orbits. This separation results from the fact 
that the local maximum of the effective potential at point 
B is greater than 0. If the angular momentum is reduced 
below = 7.58 for 5 = 2, E ^ 0.97 (Fig. W), the 
maximum B drops below and a saddle point appears. 
This saddle point corresponds to an unstable periodic 
orbit similar to the unstable periodic orbit denoted with 
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in [22, l23| , this notation is adopted in the 
present article as welL 

In the integrable case S = 1 (Schwarzschild metric) 
from such an unstable orbit "x" emanate the branches 
of the separatrix manifold, which separate the plunging 
from the non-plunging orbits. When the system is non- 
integrable, instead of the separatrix manifold, there are 
asymptotic manifolds emanating from the "x" orbit, "x" 
is formally called Lyapunov orbit (LO) [S^. Every or- 
bit crossing the border defined by a LO "escapes" , which 
in our case means that every orbit which passes the LO 
while moving towards the central object, will plunge to 
the central object. On the other hand, for certain energy 
and angular momentum the "x" orbit changes its sta- 
bility to indifferently stable and becomes the innermost 
stable circular orbit (ISCO). 

In order to locate the ISCO, it is more convenient for 
the algebraic manipulation to use the potential V — 
2 Vcs/Z, instead of the effective potential (fT3| . For 
the ISCO radius xisco the potential V and its two 
first derivatives with respect to x are equal to zero, i.e. 



dV I 

dx l^isco 



a^isco 

^^ISCO 



aS2-|a;isco 



= 0. Then we find 



3(5+ V552~~T, 



-gu\xr, 



a^isco 



* or' 

2;isco - 2d 



(14) 
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for 6 > l/\/5. For 5 < l/VS, on the other hand the ZV 
spacetime ceases to have an ISCO. From a dynamical 
point of view, below this limit we may consider the ZV 
rather as a perturbation of the Minkowski {6 = 0) than 
of the Schwarschild {S = 1) spacetime. The ISCO in the 
ZV spacetime was studied recently in [ioj , where the au- 
thors investigated the possible observational differences 
between the properties of accretion disks in the case of a 
ZV spacetime and a black hole spacetime. 



2.3. Nonlinear Hamiltonian dynamics 



By simply applying a Legendre transformation 



H = p^x^ - L 



(15) 



on the Lagrangian function we get the Hamiltonian 
function 



H = 



1 
27* 



9 Pt^Pu 



(16) 



where the momenta are given by ([7]) and p'^ = fix'^ . 
It follows immediately that the ZV system is an au- 
tonomous system (^ = ^ = 0), thus the Hamiltonian 
function is an integral of motion and equal to H = — /i/2 
(because g^^p^pi, = — ^^) 



As already mentioned in section I2.2[ the study of an 
axisymmetric and stationary system can be reduced to 
the meridian plane. Therefore, the system is reduced 
to a Hamiltonian system with two degrees of freedom. 
In the case of the Schwarzschild metric {5 = 1), spher- 
ical symmetry introduces an extra integral of motion, 
namely the total angular momentum. This integral is 
the reason why the Schwarzschild metric corresponds to 
an integrable system. In this case bounded non-plunging 
orbits oscillate in both degrees of freedom with two char- 
acteristic frequencies, denoted hereafter wi and uj^- The 
motion is restricted on a two-dimensional torus, called 
invariant torus. The type of motion on the torus de- 
pends on the ratio uji/u)2- If the ratio is a rational num- 
ber, the motion is periodic. The corresponding torus is 
called resonant and it hosts infinitely many periodic or- 
bits with the same frequency ratio. On the other hand, 
if L0i/u2 is irrational, the motion is called quasiperiodic. 
In that case, one single orbit densely covers the whole 
non-resonant torus. 

In general, resonant and non-resonant tori align around 
a central periodic orbit forming the so-called tori folia- 
tion. If we move away from the central periodic orbit 
along any radial direction, we encounter tori with vary- 
ing characteristic frequencies wi, uj2. In fact, the ratio 
ui /uj2 changes along the radial direction in the same way 
as the rational and the irrational numbers interchange 
along a real axis. 

Now, if an integrable system (like Schwarzschild) 
is perturbed, the transition from integrability to non- 
integrability follows two basic theorems, namely the 
Kolmogorov-Arnold-Moser (KAM) theorem [4l| and the 
Poincare-Birkhoff theorem [43 . 

According to the KAM theorem, for small perturba- 
tions, most non-resonant invariant tori of the integrable 
system survive deformed in the perturbed system. The 
new deformed tori are called KAM tori. 

On the other hand, according to the Poincare-Birkhoff 
theorem, from the infinitely many periodic orbits on a 
resonant torus of the integrable system only an even num- 
ber survive in the perturbed system. Half of these surviv- 
ing periodic orbits are stable and the other half unstable. 

The structure of the phase space is revealed conven- 
tionally by the use of a 2-dimensional Poincare surface of 
section. On a surface of section, the invariant tori cor- 
respond to closed curves, such curves are shown in Fig. 

(schematic) . For weakly perturbed systems the KAM 
curves on the surface of section are closed invariant curves 
around a stable fixed point in the center (like Uq in Fig. 
[5]). This structure is hereafter called the main island of 
stability. 

However, the details of the surface of section of a non- 
integrable Hamiltonian system close to resonances are 
quite different from those of integrable systems. The sta- 
ble surviving periodic orbits are depicted as stable points 
on a surface of section (filled circles in Fig. [2]) and these 
points are surrounded by other regular orbits forming the 
secondary islands of stability (Fig. [2]) . Between these is- 
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FIG. 2: A schematic representation of a surface of section of 
a non-integrable Hamiltonian system of two degrees of free- 
dom. The figure shows a Birkhoff chain at the resonance 1 /3 
between two KAM curves (homocentric cyclic curves around 
the central periodic point uq). The Birkhoff chains consists 
of 3 unstable points (open circles) and 3 stable points (filled 
circles). An island of stability is formed around each of the 
3 stable points (elliptic curves). From the unstable points 
emanate asymptotic curves (gray curves). The arrowheads 
on the asymptotic curves indicate the stable and the unsta- 
ble branch. The curved arrows indicate the flow in the phase 
space. The points w„, w„+i indicate two successive crossings 
of the quasiperiodic orbit forming the KAM curve. 6n+i is the 
angle between the position vectors R„, Rn+i (with respect 
to uo) of the points w„, w„+i respectively. 



lands lie unstable pomts (open circles in Fig. [2]), which 
correspond to unstable periodic orbits. From the unsta- 
ble periodic orbits emanate their asymptotic manifolds, 
yielding asymptotic curves on the surface of section (gray 
curves of Fig. [J). There arc two types of asymptotic 
curves: stable and unstable. A stable (unstable) asymp- 
totic manifold cannot cross itself or other stable (unsta- 
ble) asymptotic manifolds. Initial conditions on a stable 
(unstable) manifold tend asymptotically to the periodic 
orbit in the inverse (direct) flow of the time parameter. 
The latter property of asymptotic manifolds, combined 
with the non-crossing property produce on the surface 
of section oscillations of the manifolds, causing the so- 
called homoclinic "chaos" effect. Due to this, if we start 
with initial conditions within the domain crossed by man- 
ifolds, their subsequents give the impression of scattered 
points, which are a signature of chaotic motion on a sur- 
face of section. 

Besides the use of a surface of section, another tool to 
study non-integrable systems is the rotation number vg. 
The rotation number has been proved to be an efficient 
indicator of chaos [33l-l35| (see jisj for review). A sim- 
ple way to evaluate the rotation number is the following. 
First we identify the central periodic orbit (uq in Fig. 
[2) of the main island of stability. Then, we define the 
position vector 



of the n-th crossing w„ of the orbit through a surface 
of section with respect to the central periodic orbit Uq 
(dashed vectors in Fig. We then find the angle be- 
tween two successive vectors 9n+i = angle{'Rn+i,'Rn) 
(called the rotation angle) and we evaluate the rotation 
number 
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2ttN ^ 

n=l 



(18) 



In the limit N ^ oo, the rotation number corresponds 
to the frequency ratio ioxjuji. 

If we plot the rotation number as a function of the 
distance of initial conditions from the central periodic 
orbit Uq of the main island of stability along a particu- 
lar direction we obtain the so-called rotation curve. For 
integrable systems, like the Schwarzschild metric, the ro- 
tation curve is a smooth and strictly monotonic function. 
When the system is perturbed, the rotation curve ceases 
to be smooth, and it is only approximately monotonic. 
In fact, the imprints of chaos are found in the rotation 
curve, as demonstrated by several examples in section [31 



NUMERICAL RESULTS 




Uo 



(17) 



FIG. 3: The projections of geodesic orbits on the meridian 
plane ((2, p)) for 5 = 2, k = \,E = 0.95 and = 3. (a) The 
simple periodic orbit uo (red/gray curve) and the 2/3 periodic 
orbit (black curve), (b) A quasiperiodic (KAM) orbit. In 
both panels, the boundary of motion on the meridian plane 
is determined by the CZV (thick curve). 

Previous numerical investigations of the ZV system 
[30L I3TI I have led to a conjecture that this system is in- 
tegrable. This led to attempts of finding the missing 
integral of motion [3ll,[32j. In [32], however, the authors 
proved the non-existence of a broad category of integrals 
of motion for the (5 = 2 case. The parameter (5 = 2 was 
used also in [3l|, where the values i? = 0.95, ^^ = 3 and 
k = \ ((5 = M) were chosen for the numerical examples. 
The same values are applied in the subsequent section [3T] 
in order to provide a straight forward comparison with 
previous studies, while cases for different values of 8 are 
presented as well. 
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FIG. 4: (a) The surface of section 2: = (i > 0) for 5 — 2, E — 0.95 and Lz = 3. The dashed black curve determines the 
boundary of the bounded orbits. The closed curves are KAM curves of the main island of stability, while the empty phase 
space between the the main island of stability and the boundary curve is occupied by plunging orbits, (b) A magnification of 
the main island of stability. The point uo is the central fixed point. Various islands of stability forming Birkhoff chains are 
shown, labeled by their corresponding rotation numbers. 



3.1. Cases with 5 > 1 

The periodic orbits on the meridian plane are bounded 
by the CZV (section [2.2[) . Two examples of periodic or- 
bits projected on the meridian plane are given in Fig. 
[SK- The first is the simple periodic orbit Uq (red/gray 
curve in Fig. [3^), which bounces along an arc between 
two points of the CZV (thick curves in Fig. [3]). The sec- 
ond is a stable periodic orbit of the resonance 2/3 (black 
curve in Fig. [3^), which bounces also between points of 
the CZV but follows a more complex path. The path 
is even more complicated in the case of a quasiperiodic 
orbit (Fig. |3}d), because the quasiperiodic orbits cover 
densely their hosting KAM tori (section [2.3|) . 

Fig. m now, shows the surface of section 2; = (i > 0) 
corresponding to orbits like those of Fig. [3l The orbits 
shown on the surface of section plane (p, p) are in general 
non-equatorial. The bounded orbits lie inside the bound- 
ary curve p — ±1/— 2 Vcs{z = 0) (dashed curve in Fig. 
H^). Most of them are plunging orbits and only a small 
main island of stability of non-plunging orbits appears in 
the center of Fig. 0^. 

At first glance the magnified surface of section in Fig. 
|4}d seems to be quite regular, with no prominent signs 
of chaos. However, at the resonances 1/2, 3/5, 2/3 
thin islands of stability appear. These indicate the ex- 
istence of Birkhoff chains and, therefore, of chaos. In 
fact, if we zoom at the left tip of the main island of 
stability we already get a typical picture of a chaotic 
layer (Fig. [5]). In Fig. [5] the scattered points define 
chaotic regions, while the continuous curves define the 
limits of regular domains. However, the separation be- 
tween chaotic and regular regions is not so clear. For 




FIG. 5: Detail of the left tip of the main island of 
stability shown in Fig. [IJa. The arrows are showing 
the position of the islands of stability of the resonances 
11/14, 18/23, 31/40, 10/13 along the p = line. 



example, the islands of stability belonging to the reso- 
nances 11/14, 18/23, 31/40 are embedded in prominent 
chaotic layers (Fig. [5]). 

Most of the chaotic orbits of Fig. [S]are plunging orbits 
which, due to the phenomenon of stickiness |44| . remain 
for a long time close to regular non-plunging orbits, be- 
fore crossing the LO and plunge. The LO lies approxi- 
matively at the edge tip formed at the left of the main 
island of stability of Fig. [5] In fact the region between 
the main island of stability and the boundary curve (Fig. 
H^) is covered by chaotic plunging orbits. The chaotic 
orbits stick around islands of higher multiplicity, e.g. the 
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islands of stability 11/14, 18/23, 31/40 (Fig. [5|), before 
they plunge. The multiplicity is equal to the denominator 
of the prime number ratio corresponding to the rotation 
number. 




7.52 7.53 7.54 7.55 

P 



FIG. 6: The rotation number as a function of the initial con- 
ditions lying on the p = line of Fig. (5] Prominent islands 
of stability are labeled by the respective rational value of ug. 

The rotation curve along the p = line of Fig. [S] is 
shown in Fig. [6] The regions dominated by regular mo- 
tion are represented by relatively smooth segments of the 
curve, while the chaotic regions are recognized by various 
distinct parts where the rotation curve develops fluctua- 
tions. At certain smooth segments of the rotation curve, 
the rotation number is constant. These segments appear 
like "plateaus" of the rotation curve. Distinct plateaus 
correspond to islands of stability of distinct resonances 
in Fig. [5] In Fig. |5] the most prominent plateaus are 
labeled by the value of the rotation number at each re- 
spective resonance. 
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FIG. 7: The rotation curve along the line p = corresponding 
to the Fig. |4]3. The embedded figure shows a detail of the 
curve at the 2/3 resonance. 



In fact, the computation of the rotation curve is an effi- 
cient way to look for initial conditions leading to Birkhoff 
islands of stability, since the latter can be located by the 
plateaus of the rotation curve. This method is partic- 
ularly useful in the location of tiny, or narrow islands, 
such us in Fig. ^p. In this case, by scanning along 
the line p = the corresponding rotation curve (Fig. 
[S]) seems to be rather smooth (excluding the last part 
corresponding to the chaotic layer surrounding the main 
island of stability). However, a detailed scan reveals a 
number of lower multiplicity resonances (e.g. the reso- 
nances 1/2, 2/3, 3/5 shown in Fig. 2]). Lower multiplic- 
ity resonant structures (like islands of stability) are more 
prominent than higher multiplicity ones (as e.g. the res- 
onance 10/13 shown in Fig. [5]). By focusing on the lower 
multiplicity resonances the anticipated plateaus of the 
islands of stability clearly appear on the rotation curve. 
An example is shown for the plateau of the conspicuous 
2/3 island of stability in the embedded panel of Fig. [71 




FIG. 8: (a) A detail of the 2: = (i > 0) surface of section 
for (5 = 3, £ = 0.98, = 7 and fc = 1. An island of 
stability belonging to the 2/3 resonance is shown embedded 
in a chaotic layer, while the arrows indicate the position of 
the unstable periodic points of other Birkhofi' chains, (b) The 
rotation curve along the line p = corresponding to Fig. 
[8^. The plateau of the 2/3 is denoted and the arrows denote 
the behavior of the rotation curve when it crosses unstable 
periodic points. 

If we change the value of the 5 parameter to 5 = 3 
the structure of the phase space does not change signifi- 
cantly. Near the outer boundary of the main island there 
is a sticky chaotic layer of plunging orbits (Fig. [5^) and 
Birkhoff chains appear inside the main island of stability. 
The corresponding rotation curve in Fig. [SId shows large 
variations whenever it crosses initial conditions belong- 
ing to chaotic orbits. Also, the rotation curve takes the 
form of a plateau when crossing resonant islands of sta- 
bility, and it changes abruptly when crossing the unstable 
periodic points of relatively small resonances (arrows in 
Fig. [5}d). Thus, in the case of ^ = 3 as well both detect- 
ing methods, i.e. the surface of section and the rotation 
number, indicate that the ZV metric is not integrable. 

If we further increase the value of 5, the overall picture 
of the phase space is qualitatively similar to the cases 
already presented. However, what changes is the posi- 
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FIG. 9: (a) A detail of the z — {z > 0) surface of section for 
S = 4, E = 0.95, Lz = 3 and fc = 1. (b) The rotation curve 
along the line p = corresponding to Fig. [9ji. 



tion of the unstable point "x" (section [2. 2 1) , which moves 
further away from the central anomaly p — z — 0. This 
is expected, since the position of the ISCO moves away 
as well (first equation of eqs. (|14p ). Thus, for increas- 
ing S the main island of stability moves all together to 
larger distances p. For example. Fig. IH] shows a detail 
of the phase space for 5 = 4, i? = 0.95, — 3, where 
the whole island structure has been shifted further away 
from p = z = with respect to Figs. [SJ [HI [E 




4.805 4.81 4.815 4.82 4.825 4.8.1 4.8,15 




4.805 4.81 4.815 4.82 4.825 4.83 4.835 



FIG. 10: (a) A detail of the 2: = (i > 0) surface of section 
for S = 3/2, E = 0.98, = 5.7 and = 1. (b) The rotation 
curve along the line p — corresponding to Fig. llOb . 

In all previous cases the value of S was a natural num- 
ber. However, the structure of the phase space seems not 
to change even if (5 is a positive real (not integer) num- 
ber. For instance, if (5 = 3/2, then we find again a chaotic 
layer around the main island of stability (Fig. [TU)) . and 
all the complexity seen in the previous examples as well. 
In particular, islands of stability are either embedded in 
a prominent chaotic layer or enveloped between KAM 
curves (Fig. [TOk). The resonances corresponding to these 
islands are found in Fig. [TUb by the use of the rotation 
number. Fig. [TUb is a good example to note again that 
lower multiplicity islands of stability are more prominent 
than the higher multiplicity ones. This fact is the reason 
why we expect [Til . ITs | that the lower multiplicity islands 
of stability are good candidates for detecting non-Kerr 



compact objects by the analysis of gravitational waves 
coming from EMRIs, even if the inspiraling smaller com- 
pact object might cross infinite resonances in a bumpy 
black hole spacetime background during its inspiral. 



3.2. Cases with 5 <1 





0.814 0.816 0.818 0.S2 0.822 



FIG. 11: (a) A detail of the 2: = (i > 0) surface of section 
for S = 1/2, E ^ 0.99, = 1.8 and k = 1. (b) The rotation 
curve along the line p = corresponding to Fig. Illb . 

When (5 < 1, the ZV spacetime corresponds to a more 
prolate central object than the respective Schwarzschild 
black hole. Nevertheless, for the prolate case ^ = 1/2 
the structure of the phase space is similar to the oblate 
cases examined in section [3.11 There is a main island of 
stability, and around it a chaotic sea of plunging orbits. 
In particular, the detail of the surface of section z = 
(i > 0) (Fig. [TTk ) presents a similar phase space struc- 
ture near the sticky chaotic layer around the main island 
of stability as in the oblate cases shown in section 13.11 
Thus, we find again the 3/4 Birkhoff chain, from which 
we can see clearly one of the respective islands of stability 
embedded in the chaotic sea. Islands of higher multiplic- 
ity lie near the border of the chaotic layer and of the 
main island of stability. Thus, it seems that the ZV sys- 
tem has a similar dynamical behavior independently of 
whether the corresponding central object is more oblate 
or prolate than a Schwarzschild black hole. 

However, if we reduce further the parameter 6, for in- 
stance to S = 1/3, then the dynamical behavior changes 
significantly. This change is expected for S < 1/ \/5, be- 
cause as noted in the last paragraph of section 12.21 the 
ZV spacetime loses a property, namely the existence of 
an ISCO, which is a characteristic of a black hole space- 
time. In Fig. [1^ the surface of section shows a main 
island of stability surrounded by a chaotic sea. Due to 
the fact that the CZV is a closed curve (bottom right 
embedded panel in Fig. [T2|) . the orbits of the chaotic 
sea cannot plunge towards the central anomaly. As a 
result, even if the chaotic sea is prominent in size, no 
strong chaotic regions can be seen inside the main island 
of stability. This fact is confirmed both by the surface of 
section and by the rotation curve (main panel and right 



9 




FIG. 12: The z = {i > 0) surface of section for 5 = 1/3, 
E = 0.99, I/z = 0.4 and k — 1. In the lower right embedded 
figure is shown the corresponding CZV. In the upper right 
embedded figure is shown the rotation curve along the p — 
for the segment corresponding to the main island of stability. 



up panel in Fig. [12] respectively) . Thus, the structure 
of the main island of stability seems not to change sig- 
nificantly for S < (compare the surface of section 
and the rotation curve in Fig. [T2]with that in Figs. HIE] 
respectively). 



4. CONCLUSIONS AND DISCUSSION 

It has been shown by clear numerical examples that 
the ZV metric corresponds in general to a non-integrable 
system, contrary to the results of previous works [sol . 
l3l| . The only known ZV spacetimes that correspond to 
integrable systems, are the Minkowski spacetime {6 = 0) 
and the Schwarzschild spacetime (i5 = 1). In order to 
find imprints of the ZV non-integrability the methods of 
the surface of section and of the rotation number were 
employed. 

In fact, in cases like the ZV system, we can only guess 
where to zoom on a surface of section to discover chaotic 
layers or other imprints of chaos. However, in such 
cases of nearly integrable Hamiltonian systems, chaos is 
strongly correlated with the remnants of the destroyed 
resonance tori (section |2.3|) . Therefore, a tool, like the 
rotation number, which detects resonances is very useful. 

Moreover, the rotation number is a reliable indicator of 
chaos. Thus, even if we do not use a surface of section or 
another method for detecting chaos, the rotation curves 
themselves are sufficient to show that chaos exists. This 
is interesting also from an observational point of view, 
because the rotation number is the ratio of the two main 
frequencies of a non-plunging geodesic orbit. Therefore, 
the rotation number is an adequate tool for detecting 
phenomena associated with chaos in gravitational wave 
signals. This has been already studied in [Ulli^, where 



the rotation number was used in order to detect devia- 
tions from the Kerr metric in a case of an EMRI into a 
bumpy black hole spacetime background. 

It has been shown that the ZV spacetime has a similar 
dynamical behavior in all cases, i.e. independently of 
whether the spacetime corresponds to a more oblate or 
prolate central compact object, as long a.s S > l/Vb. 
In particular, all the orbits belonging to the chaotic sea 
which surrounds the main island of stability are plunging. 
When, < 5 < 1/V5 then the ISCO disappears and 
the orbits of the chaotic sea cease to be plunging. Thus, 
6 = 1/ ^/E is a marginal value, where the ZV metric family 
changes its behavior. 

Finally, the study of bumpy black holes ISCOs has an 
interesting aspect regarding spin measurements of Kerr 
black holes. If we suppose that ZV is describing the 
spacetime around a compact object, then in the inter- 
val 1/^/5 < S < 1 we get ISCO positions that correspond 
to the Kerr black hole ones. This fact implies that if we 
evaluate the central compact object spin by calculating 
the ISCO position like in [45| without having tested the 
Kerr hypothesis first, then these measurements might be 
misleading. 
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APPENDIX A: NUMERICAL ACCURACY 
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FIG. 13: (a) The relative error 



0.5 1. 1.5 2. 2.5 3. 



of the Lagrangian 



function L, which is a constant, as evaluated at the A'^th sur- 



face of section, (b) The relative error | '^'^'^ \ between the 

Lagrangian function two consecutive integration steps at the 
A'^th surface of section. The red (gray) points corresponds to 
a chaotic orbit with initial condition p = z = 0, p = 7.518, 
while the black to a regular orbit with initial condition p = 
z = 0, p = 7.548. Both initial conditions lie on the surface of 
section shown in Fig. [Sj 
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For integrating the geodesic orbits in the ZV system 
the 5th order Cash-Karp Runge-Kutta was applied with 
a step size controller developed by the author. The con- 
troller's function is to keep the relative error | — ^-^f^— | 
of the evaluated Lagrangian Lt, Lr+dr functions between 
two consecutive integration steps r, t + dr below a tol- 
erance value. Even if such a control can secure the accu- 
racy between two consecutive integration steps, it cannot 
guarantee the accuracy of the overall calculation. The in- 
tegration accuracy is checked by the relative error | ^^^"^ | , 
where the evaluated Lagrangian function L^- is compared 
with the theoretical one L. 



The overall time evolution of the relative error | | 
at the Nth surface of section crossing is shown in Fig. 
[T3k for a regular orbit (black dots) and for a chaotic 
orbit (red/gray points). The evolution of the relative 
error | ^^^^ — - \ for the same orbits is shown in Fig. [13b . 
Although tlie relative error at consecutive steps appears 
rather small, the overall error shows that both chaotic 
and regular orbits exhibit slow numerical drift. However, 
in both cases this drift is very small compared to the 
scale of studied phenomena. 
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